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Abstract 



Numerical simulations of convection driven rotating spherical shell dynamos have often been performed with rigid boundary con- 
y— I [ditions, as is appropriate for the metallic cores of terrestrial planets. Free-slip boundaries are more appropriate for dynamos in 
other astrophysical objects, such as gas-giants or stars. Using a set of 57 direct numerical simulations, we investigate the effect of 
free-slip boundary conditions on the scaling properties of heat flow, flow velocity and magnetic field strength and compare it with 
earlier results for rigid boundaries. We find that the nature of the mechanical boundary condition has only a minor influence on 
the scaling laws. We also find that although dipolar and multipolar dynamos exhibit approximately the same scaling exponents, 
there is an offset in the scaling pre-factors for velocity and magnetic field strength. We argue that the offset can be attributed to the 
differences in the zonal flow contribution between dipolar and multipolar dynamos. 
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p ^ 1. Introduction 

Numerical simulations of dynamos in geometries appropri- 
C"| .ate for the cores of terrestrial planets have greatly enhanced our 
Q-funderstanding of the complex magnetic field behavior observed 
q .in these objects, with possible implications for a broader class 
!— j 'of dynamos in astrophysical objects (Jones, 201 1). One of the 
c/3 .major drawbacks of such simulations is that there is order of 
^3 -magnitude disagreement between the natural and the numeri- 
cally accessible values of several control parameters. For in- 
t-H 'stance, in numerical simulations of the geodynamo, the Ekman 
number - a nondimensional measure of the importance of vis- 



cous effects as compared to the Coriolis effects - is usually five 



to ten orders of magnitude larger than the expected realistic val- 

"ues. 

* . One way to tackle this disparity is to infer asymptotic scal- 
'ing laws from a sufficient number of numerical results. Such 
^NJ numerical scaling laws can then be extrapolated to realistic 
^ — I parameter regimes and compared with the observational data. 
Christensen and Aubert (2006) (hereafter "CA6") used a bat- 
tery of numerical simulation results to derive scaling rela- 
tions for heat transfer, convective velocity, and magnetic field 
strength. Their scaling relations hold over several orders of 
magnitude of the relevant control parameter. Using these scal- 
ing relations, CA6 predicted magnetic field strengths inside the 
Earth's and Jupiter's core and found reasonable agreement with 
observationally constrained values. Later, Christensen et al. 
(2009) showed that these scaling laws are also in good agree- 
ment with magnetic fields observed in fast rotating low-mass 
stars. Takahashi et al. (2008) and Aubert et al. (2009) inde- 
pendently re-enforced the scaling laws put forth by CA6. 
Olson and Christensen (2006) derived scaling laws specifically 
for the dipole moment of planetary dynamos which show an 
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order-of-magnitude agreement with the observed dipole mo- 
ments of solar system planets. Christensen (2010) reviews ear- 
lier scaling laws for planetary magnetic field based on heuristic 
arguments and compares them with the numerically established 
scaling relations. 

The mechanical boundary conditions may play an im- 
portant role in the dynamo mechanism. Dynamos which 
operate in planets with a solid mantle are usually mod- 
eled with rigid boundaries (Kageyama and Sato, 1995; 
Glatzmaier and Roberts, 1995a,b). While a true free surface 
is difficult to model, a free-slip condition (i.e. assuming zero 
shear stress at an undeformable spherical boundary) is a much 
better approximation than a no-slip condition for the surface 
of gas- and ice-giant planets or stellar convection zones. Rigid 
boundaries are associated with viscous (Ekman) boundary 
layers, which have a damping influence on the development 
of strong axisymmetric flows that are found in free surface 
flows. Kuang and Bloxham (1997, 1999) argued that even for 
the geodynamo a free-slip condition may be a better choice, 
because the Ekman layers in the models are much thicker than 
the very thin layers in the Earth's core. Aubert et al. (2001) 
demonstrated in rotating liquid Gallium experiments that at 
small Prandtl numbers - the ratio of kinematic viscosity to 
thermal diffusivity - flows with rigid boundaries show features 
similar to those with free-slip boundaries, i.e. dominant zonal 
flows. Miyagoshi et al. (2010) also found strong zonal flows in 
low Ekman number 10~ 7 ) rigid boundary geodynamo sim- 
ulations. Hence, studies of dynamos with free-slip boundaries 
have wide ranging applications. Following these arguments, 
researchers have modeled Uranus' and Neptune's multipolar 
magnetic field (Stanley and Bloxham, 2004, 2006), Mercury's 
weak magnetic field (Stanley et al., 2005), and Saturn's un- 
usually axisymmetric magnetic field (Stanley, 2010) using 
free-slip boundaries. 

Dynamos with free-slip boundaries also exhibit bista- 
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bility: the morphology of the dynamo generated mag- 
netic field depends on the initial magnetic field configu- 
ration (Simitev and Busse, 2009, 2012; Sasaki et al., 2011; 
Schrinner et al., 2012; Gastine et al., 2012). Grote and cowork- 
ers (Groteetal., 1999; Grote and Busse, 2000; Grote etal., 
2000) have employed free-slip boundaries in their dynamo 
models, and found a wider spectrum of magnetic field geome- 
tries than what has been reported for rigid boundaries. 

A direct comparison of the effects of different mechanical 
boundary conditions on the dynamo has rarely been made. 
Christensen et al. (1999) reported results for both kinds of 
boundary conditions for a limited number of cases and found 
that the large scale magnetic field is similar for both cases. 
Recently, Schrinner et al. (2012) have analyzed many dynamo 
simulations with rigid, free-slip and mixed (rigid at inner and 
free-slip at outer boundary) boundary conditions and report a 
difference in magnetic field amplitude of dipolar and multipo- 
lar dynamos. Following this study, we specifically focus here 
on deriving scaling properties for heat transport, velocity, and 
magnetic field strength. We compare our findings with earlier 
rigid boundary systems for which extensive modeling results 
are available in a broad range of control parameters (CA6). This 
exercise helps us in isolating the effect of mechanical boundary 
condition. 

2. Dynamo model 

2.1. MHD equations 

Our numerical set-up consists of a spherical shell which ro- 
tates along the z-axis and which is bounded by inner radius r,- 
and outer radius r„. The aspect ratio r,/r is 0.35. A linear vari- 
ation of gravity with radius is assumed. Following previous pa- 
rameter studies, we non-dimensionalize the magnetohydrody- 
namic (MHD) equations by using the shell thickness r - r, = D 
as the reference length scale and D 2 /v, where v is the fluid 
viscosity, as the time unit. The magnetic field B is scaled by 
-y/p///lf2, where p is the constant fluid density, ju is the magnetic 
permeability, A is the magnetic diffusivity, and Q. is the rotation 
rate. The temporal evolution of velocity u, temperature T, and 
magnetic field B is governed by the MHD equations under the 
Boussinesq approximation 

E |^ + u • Vu - V 2 uj + 2z x u + VP = 

Ra — + — (VxB) xB, (1) 

r„ Pm 

^ +u . V r = i-v 2 r, (2) 

at Pr 

y=Vx(uxB) + iv 2 B, (3) 
at Pm 

V ■ u = 0, (4) 

V ■ B = 0. (5) 

This system of equations is governed by several nondimen- 
sional control-parameters: Ekman number E = v/Qd 2 ; the 
modified Rayleigh number Ra = ag ATD/vQ., where g a is 



gravity at the outer boundary and a is the thermal expansivity; 
magnetic Prandtl number Pm = v/A\ Prandtl number Pr = v/k, 
k being the thermal conductivity. Ra is related to the conven- 
tional Rayleigh number Ra through Ra = RaE/Pr 

We assume free-slip mechanical boundaries at both inner and 
outer radius. The magnetic field matches a potential field at 
both boundaries. A fixed temperature contrast AT is maintained 
between the top and the bottom. 

2.2. Numerical method 

Equations (1-5) are numerically solved using the MagIC 
code (Wicht, 2002). Velocity and magnetic field are first sepa- 
rated into toroidal and poloidal components as 

u = V x u T f + V x V x u P f, 
B = Vxfij-r + Vx V xBpf. 

The scalar potentials ujp and Bjp, along with temperature T 
and pressure P, are further expanded using spherical harmonics 
in the 6 and <p directions and the Chebyshev polynomials in the 
radial direction. N r and l max are the maximum degree of the 
Chebyshev polynomials and the spherical harmonic functions 
used in this expansion. For all the simulations considered here, 
41 < N r < 73 and 64 < l max < 170. The simulations are 
run for at least one magnetic diffusion time (D 2 /if) to ensure a 
statistically stationary state. 

2.3. Diagnostic parameters 

We employ several diagnostic parameters to quantify the 
model results in terms of spatially and temporally averaged val- 
ues. The Rossby number Ro quantifies the mean velocity u and 
is defined by Ro = u/(DQ). Following CA6, we also intro- 
duce the local Rossby number Ro/ which is a more appropriate 
measure than Ro to characterize the ratio of the inertial and the 
Coriolis forces. A typical flow length scale can be calculated 
based on the mean spherical harmonic degree / 

L = ^ l (»t ■ u/>, 

where (...) denote time average and u/ is velocity component at 
degree I. Roi is then defined as Rot - Ron/l,,. 

The mean magnetic field strength is often expressed by the 
Elsasser number A = B 2 m J ^jppAQ.. However, for ease of com- 
parison with CA6, we use the Lorentz number Lo = y/EA/Pm 
to quantify the magnetic field strength. The field geometry at 
the outer boundary surface is characterized by its dipolarity fa p . 
It is the ratio of the magnetic energy of the dipole to the total 
magnetic energy at the outer boundary surface. 

The Nusselt number Nu is a ratio of total heat transported 
from the inner shell to the outer shell to the conducted heat. It is 
expected that different diffusivities play a minimal role in deter- 
mining the large scale properties of the dynamo systems. This 
motivated CA6 to use the 'diffusionless' parameters Ro and Lo 
and also to define a modified Nusselt number Nu* which does 
not involve k. Nu* is related to the conventional Nusselt number 
Nu via Nu* = (Nu - \)E/Pr. In addition, the heat flux from sur- 
faces of astrophysical objects is a much more meaningful and 
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accessible quantity than the temperature difference between the 
inner and outer boundary of the convection zone. CA6 defined 
a heat flux based Rayleigh number Ra* Q which incorporates the 
advected heat flux rather than the temperature contrast. Ra* Q is 
related to Ra through Ra* Q = Ra(Nu - l)E 3 /Pr 2 . 

3. Results 

We have built up a data set of 57 dynamo simulations with 
free-slip boundaries: 40 cases by us (see Table A. 2) and 17 
cases adopted from Schrinner et al. (2012). All of the cases in 
this data set have the same physical set-up as described above. 
Also, following CA6, we report and analyze only those dynamo 
simulations which have Nu > 2 to ensure a vigorous enough 
convection that fills the full volume of the spherical shell. 

3.1. Bistability 

Bistability is a phenomenon in which a system shows two 
different dynamo solutions, i.e. dipolar and multipolar, for 
the same set of control parameters, but with different initial 
conditions for the magnetic field. For dynamos with free- 
slip boundaries, Simitev and Busse (2009) found two distinct 
dipolar and multipolar dynamo branches, at least in some pa- 
rameter range. Recently, Sasaki et al. (2011), Schrinner et al. 
(2012), and Gastine et al. (2012) have also observed bistabil- 
ity in their dynamo simulations with free-slip boundaries. Such 
bistable branches have rarely been observed in dynamos with 
rigid boundaries, e.g. CA6 reported a single case in which they 
found bistable states. 

We plot in Fig. 1 the dipolarity f^ p at the outer surface 
(r = r a ) versus the local Rossby number Roi (for our simula- 
tions). In the dipolar branch (filled data points) the magnetic 
field is dominated by the dipole {Ja p > 0.4) and in the mul- 
tipolar branch (empty data points) the dipole is much weaker 
ifdip < 0.05). The dipole branch is limited to cases with 
Roi < 0.2. However, the multipolar branch exists for a broad 
range of Rot. The highly supercritical dipolar case at Roi * 0.2 
(E = 1 x 10~ 4 , Pm = 0.5 and Ra = 4000) was run for two 
magnetic diffusion times without an indication of a dipole col- 
lapse, although we can not exclude that the field could change 
to multipolar in long run. Dipolar dynamo at such high Ro/ has 
not been reported yet. Simulations that settle down to differ- 
ent dynamo states depending on the initial magnetic field are 
marked with a "+" in Fig. 1 (difficult to discern on the multipo- 
lar branch due of clustering). Note that other dipolar dynamos 
could show bistability but we did not explore all of our dipolar 
cases for such behavior. 

For rigid boundary dynamos, CA6 reported that multipolar 
dynamo solutions are not observed for Ro/ < 0.1. But results 
from Schrinner et al. (2012) and our findings suggest that this 
is not the case for free-slip boundaries. We also found few dy- 
namo solutions, which have Rot < 0.1, but settle to a multipolar 
solution despite having initial dipolar magnetic field. For exam- 
ple, we only found a multipolar solution for E — 3 x 10~ 4 and 
Pm < 1.5. This demonstrates that, depending on the control 
parameters, only the multipolar dynamo branch can be stable 
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Figure 1 : Dipolarity at the outer boundary versus the local Rossby number. The 
data points carrying "+" marker are bistable states. Filled (empty) symbols are 
dipolar (multipolar) dynamos. The symbol shapes represent the Ekman number 
and the corresponding value is given in the upper-right-corner box. 



in some situations. This is in agreement with earlier results 
(Simitev and Busse, 2009, 2012) which showed that bistability 
is in-fact a function of P, Pm and E. They used volumetric 
heating while we use a fixed temperature contrast to drive con- 
vection. The bistable behavior we observed generalizes their 
findings. 

One of the characteristic features of rotating spherical shell 
convection with free-slip boundaries is the development of 
strong axisymmetric zonal flows. In rigid boundary systems, 
the Reynolds stresses, which arise due to a statistical correlation 
between the radial and azimuthal flow component (in cylindri- 
cal co-ordinates), are balanced by the bulk viscosity and the 
Ekman layer friction near the outer boundary. In case of free- 
slip boundaries, the Ekman layers are absent and zonal flows 
can thus saturate at much higher vigour. In dynamo models, 
Maxwells stresses also affect the zonal flows; these stresses are 
potentially higher in dipolar dynamos that have higher mag- 
netic field strength than multipolar ones at the same control pa- 
rameter values (Browning, 2008). One argument for the essen- 
tial role of zonal flows for bistability is that an initial dipolar 
magnetic field inhibits the growth of zonal flows via Maxwell 
stresses. In case of a multipolar initial condition, strong zonal 
flows can develop, which in turn suppress the development of 
dipolar magnetic fields. This mechanism allows multipolar 
magnetic fields even for Ro/ smaller than 0.1. The resulting 
zonal flow structure of a bistable state is shown in Fig. 2. It por- 
trays a weak thermal-wind driven zonal flow in the dipolar case 
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Figure 2: Non dimensional radial magnetic field at outer boundary and zonal 
flow (azimuthally averaged) in a meridional section for a dipolar (top row) and 
a multipolar (bottom row) dynamo. These states are obtained at E = 3 X 10~ 5 , 
Pm = 0.5, and Ra = 900. The radial magnetic field is truncated at 50% of the 
maximum in order to highlight the magnetic field structures. 



and a nearly three times stronger geostrophic zonal flow in the 
multipolar case. The non-geostrophy of the flow in case of the 
dipolar dynamo emphasizes that the Lorentz force is actively 
participating in the force balance of the dynamo mechanism, 
while this is not the case in the multipolar dynamo. As also 
pointed out before by Schrinner et al. (2012). 

3.2. Nusselt number scaling 

Figure 3 shows that the modified Nusselt number Nu* scales 
very well with the flux-based Rayleigh number Ra* Q in the same 
way irrespective of whether the dynamo is dipolar (filled sym- 
bols) or multipolar (empty symbols). A best-fit line to this data 
set reveals a relation Nu* - 0.061 Rcl*q 52 , which agrees well 

with the scaling Nu* - 0.076 Ra* Q - 53 (dashed line) found by 
CA6 for dipolar dynamos with rigid boundaries. For convection 
in rotating shells without magnetic field, Christensen (2002) 
suggested a possible asymptotic scaling Nu* = 0.077 Ra*^ 9 
in the limiting case of E — > 0. All these scaling relations are 
very close to each other, which suggests that neither the mag- 
netic field nor the mechanical boundary conditions have sub- 
stantial effect on the scaling behavior of heat transport in ro- 
tating spherical shell convection. However, we note that for 
larger values of the Rossby number than the ones considered 
here, in a regime where inertia dominates over the Coriolis 
force, the power-law scaling between Nu* and Ra* Q breaks down 
(Schmitz and Tilgner, 2010; King et al., 2010). 

3.3. Rossby number scaling 

In Fig. 4 we plot Ro as a function of Ra* Q . The data points 
are slightly more scattered as compared to Fig. 3, but a con- 
sistent scaling is nonetheless evident. Moreover, a somewhat 
different scaling for dipolar and multipolar dynamos is visible 
as demonstrated by the two different solid lines. These lines are 
Ro = 0.73 Ra* Q 039 (dipolar) and Ro = 1.78 Ra* Q 0M (multipo- 
lar). The scaling reported by CA6 for dipolar rigid boundary 
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Figure 3: Modified Nusselt number versus heat flux based Rayleigh number. 
The data symbols carrying a thick black dot are adopted from Schrinner et al. 
(2012). The solid-line is a line-fit to data and the dashed-line is the scaling 
reported by CA6 for dipolar rigid boundary dynamos. 



dynamos is Ro = 0.85 Ra* Q (dashed-line), which agrees with 
our dipolar dynamo scaling. 

Similar to CA6, the scatter in Fig. 4 can be reduced 
to some extent by assuming an additional Pm dependence. 
A two-parameter least-square-optimized fit provides Ro = 
0.99Ra* Q 0A1 Pm' 01 (dipolar) and Ro = 2.44 Ra* Q 0A1 PnT 014 
(multipolar). This optimization reduces the standard error by 
48% in the dipolar scaling and 20% in the multipolar scaling 
(see Table 1). We also considered the Ekman number as addi- 
tional parameter for improving the fit, but, as observed by CA6, 
the resulting exponents are rather small as compared to Ro and 
Pm exponents. Hence, we discard a dependence on E in our 
scaling analysis. The Pm exponent is small and appears to de- 
pend on the nature of the magnetic field. The latter could be 
an artifact of relatively small size of the data set, especially for 
dipolar dynamos. In fact, similar to CA6, assuming a scaling 
of the form Ro oc (Ra*QPm^ 1 ^) , a good fit is obtained for both 
dipolar and multipolar cases, with somewhat different values 
for a (see Table 1). This simplified form is shown in Fig. 5. 

As described in Sec. 3.1, zonal flows are stronger in ro- 
tating convective shells with free-slip mechanical boundaries. 
In the corresponding MHD systems, dynamos with multipo- 
lar magnetic fields will have stronger zonal flows as compared 
to those with dipolar magnetic fields. This effect is visible 
in Fig. 6 whose top panel shows Ro non - zona i versus Ra* Q , and 
the bottom panel shows Ro zona t versus Ra* Q . Ro non - zona i and 
Ro zona i are calculated by considering the rms velocity exclud- 
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Figure 4: Rossby number versus the heat flux based Rayleigh number. The 
two solid lines are best-fitting lines to dipolar and multipolar dynamos. The 
dashed-line represent the scaling reported by CA6 for dipolar rigid boundary 
dynamos. 
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Figure 5: A simplified form of Rossby number scaling incorporating a Pm 
dependence. 
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Figure 6: Non-zonal Rossby number in (a) and zonal Rossby number in (b) 
versus the heat flux based Rayleigh number. 



ing the axisymmetric zonal-flow component and the rms ve- 
locity of only the axisymmetric zonal flow, respectively. The 
scaling in Fig. 6a is Ro non - zom i = 1.37 Rci*q )M and in Fig. 6b is 
Ro ZO nai = 0.29 Ra* Q QM (dipolar) and Ro ZOIM , = 0.73 Ra* Q QA (mul- 
tipolar). The zonal Rossby number of dipolar dynamos is con- 
sistently lower than that of corresponding multipolar dynamos, 
while the non-zonal Rossby number is almost unaffected by the 
field topology. This difference in zonal flows explains the off- 
set in the Rossby number scaling in dipolar and multipolar dy- 
namos seen in Fig. 4 and 5. Note that as compared to the dipolar 
branch the scatter in the multipolar branch of Fig. 6b is large. 
This could be due to the fact that unlike the dipolar branch the 
multipolar branch is a blend of dynamos which have qudrupo- 
lar, octopolar, and sometimes even higher order modes as the 
most dominating magnetic mode. Since the Maxwells stresses 
are dependent on the magnetic field geometry, the zonal flows 
saturate at many different levels. 

3.4. Magnetic field scaling 

The Ohmic dissipation time T mag , which is the ratio of mag- 
netic energy and Ohmic dissipation, is a function of the typical 
length scale of the magnetic field. As the magnetic Reynolds 
number is increased, the magnetic field becomes smaller scaled, 
and, since small scales are associated with faster time scales, the 
Ohmic dissipation time scale decreases. This qualitative argu- 
ment was verified by Christensen and Tilgner (2004) in rigid 
boundary spherical shell dynamos. They showed that T mag 
(normalized with magnetic diffusion time) is approximately 
inversely proportional to the magnetic Reynolds number Rm. 
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Figure 7: Ohmic dissipation time versus the Rossby number. Solid line repre- 
sents a best fit line, while the dashed line represents T' mag oc 1 /Ro. The inset- 
figure contains bistable pairs connected by solid lines. 




Figure 8: Corrected Lorentz number versus the heat flux based Rayleigh num- 
ber. The two solid lines are best-fitting lines to dipolar and multipolar dynamos. 
The dashed-line represent the scaling reported by CA6 for dipolar rigid bound- 
ary dynamos. 



When T' mag is the Ohmic dissipation time expressed in units of 
rotation period of the spherical shell, this inverse relation trans- 
lates to r' mag oc 1 /Ro. In Fig. 7 we plot -f versus Ro. A best-fit 
line to this data set suggests r' mag oc 1 /Ro 0& . Since the scatter in 
Fig. 7 is substantial, the difference between the exponents -0.8 
and -1 may not be very significant. Christensen (2010) have dis- 
cussed a more complex scaling for T mag and report a marginal 
improvement in the quality of the fit. Although the inset fig- 
ure shows a small decrease in r' mag for bistable states when the 
magnetic field is multipolar, the scalings for dipolar and multi- 
polar dynamos appears to follow the same trend. Moreover, if 
we plot t' versus Ro zona i or Ro non - zona i (not shown), then the 
scatter is increased as compared to Fig. 7. It highlights that the 
important parameter in the context of ohmic dissipation is the 
total Rossby number, which incorporates the zonal-flow contri- 
bution. 

CA6 argue that the magnetic field strength might be deter- 
mined by the power available to balance the Ohmic dissipa- 
tion. Following this argument, the Lorentz number should be 
accordingly corrected by the Ohmic fraction /„/„„ which is the 
ratio of Ohmic dissipation and the power generated via buoy- 
ancy forces. In Fig. 8, we plot the corrected Lorentz num- 
ber versus the flux-based Rayleigh number Ra* Q , which is a 
non-dimensional measure for the power generated by the ac- 
tion of buoyancy forces (CA6). A best fit is obtained by 
LolflJl = 1.09tf fl * G - 37 (dipolar) and Lo/f^ = 0.65 Ra* Q 035 
(multipolar). The dipolar scaling in Fig. 8 is close to the 
rigid boundary dipolar scaling Lo/f 1 ^ = 0.92 Ra* Q 034 (dashed 



line) reported by CA6. Furthermore, a two-parameter opti- 
mized fit for dipolar dynamos is Lo/f^ = 0.72 Ra* Q 033 Pm 0M 
and for multipolar dynamos is Lo/f^ = 0.51 Ra* Q oii Pm on . 
The inclusion of Pm reduces the standard error by 67% (dipo- 
lar) and 39% (multipolar). Again, assuming a simplified form 
L°/folm 00 (P^QP'^^y 3 ^ the quality of the fit is hardly reduced 
(Table 1, Fig. 9). A cursory inspection of Figs. 8 and 9 suggests 
that the dipolar and multipolar scalings are almost the same ex- 
cept for an offset in the pre-factor by » 8/5. 

An inverse relation of r' and Ro translates to Lo I f l [ 2 oc 

mag ' J ohm 

^Roq/Ro. If we now substitute the Ro scaling from Fig. 4 in 

the previous relation, then an offset of as 8/5 is indeed expected 
for dipolar and multipolar scaling. The t' oc 1 /Ro argument 
therefore supports the offset in the scaling observed in Fig. 8 to 
a good extent. Schrinner et al. (2012) have also reported similar 
shift in the Lo/fL, scaling. They, however, qualitatively argue 
that the offset in the scaling is due to decrease in the / /„„ in mul- 
tipolar dynamo cases. More analysis is required to conclusively 
demonstrate the reason for such an offset. 

4. Discussion and conclusions 

In this article we investigated the effect of free-slip mechan- 
ical boundaries on various scaling laws in spherical shell dy- 
namos. We compared the inferred scaling laws with earlier re- 
ported scalings for rigid boundary dynamos. 
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Figure 9: A simplified form of corrected Lorentz number scaling incorporating 
a Pm dependence. 



Table 1 : The various scaling laws inferred 
correlation coefficient R and the Standard 
square-root of number of data points). 



from our study, along with the cross- 
error (standard deviation divided by 
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0.9916 
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0.9553 
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0.9954 


0.0078 
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0.9952 


0.0078 




0.9941 


0.0064 




0.9975 


0.0039 


Lo 1 fl'l =0. 52 {Ra* Q Pm l I 3 f 33 


0.9975 
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We observed Instability, i.e. dipolar and multipolar dynamos 
coexisting for same control parameters. This agrees with the 
earlier findings (Simitev and Busse, 2009, 2012; Sasaki et al., 
2011; Schrinner et al., 2012; Gastine et al., 2012) and rein- 
forces the importance of free-slip boundaries and zonal flows 
for this phenomenon. 

The modified Nusselt number scales as Nu* - 0.061 Ra* Q 052 
which is in very close agreement with the scaling for rigid 
boundary dynamos (Christensen and Aubert, 2006) and non- 
magnetic convection in spherical shell with free-slip bound- 
aries (Christensen, 2002). This highlights that boundary con- 
ditions or magnetic field do not effect the scaling behavior of 
heat-transport in convective spherical shells. 

The Rossby number scales as Ro = 0.73 Ra* G ° 39 for dipolar 
and as Ro = 1.78 Ra*Q M for multipolar dynamos. The offset in 
the scaling of dipolar and multipolar dynamos can be attributed 
to different zonal flow characteristics: zonal flow quenching 
is stronger in a dipolar magnetic field configuration as com- 
pared to a multipolar configuration. This results in a smaller 
pre-factor in the Rossby number scaling for dipolar dynamos. 

1 /2 

The corrected Lorentz number scales as Lo f ! = 

' ■' ohm 

1.09 Ra* Q 031 for dipolar dynamos and Lo/f^ = 0.65 Ra* Q 035 
for multipolar dynamos. The exponents are almost identical but 
the pre-factors differ. We investigated the origin of such shifted 
scaling and found that, using the scaling for Rossby number 
(inferred from our data-set) and an inverse relationship between 
ohmic dissipation time and magnetic Reynolds number, parallel 

1 /2 

and shifted scalings for Lo/fL, are indeed expected. The ob- 

1 II 

served and the expected ratio of dipolar and multipolar Lolf! 
scalings were in good agreement. 

Similar to Christensen and Tilgner (2004) and 
Christensen and Aubert (2006), we also observed that a 
small dependence on the magnetic Prandtl number improves 
the quality of the scalings, especially in the dipolar dynamo 
cases. However, results form the Karlsruhe laboratory dy- 
namo experiment motivated Christensen and Tilgner (2004) 
to conjecture that such small Pm dependence might disappear 
when Pm « 1. The lowest Pm in our study is of order unity, 
which makes it difficult to ascertain this conjecture. So far, our 
free-slip simulations and the rigid boundary simulation results 
of Christensen and Aubert (2006) support scalings which have 
some Pm contribution. Simulations which attain Pm « 1 will 
shed more light on this issue. 

Christensen et al. (2009) showed that the scaling derived 
from Boussinesq dynamo simulations with rigid boundaries can 
predict the magnetic field strength in rapidly rotating stars. Ob- 
jects such as stars and giant planets have free-surface flows, 
very high density stratification, and probably fully convective 
interiors. Similar scaling of physical properties despite such 
drastic physical differences is puzzling. Our analysis shows that 
the boundary conditions do not substantially effect the scaling 
behavior of the rms velocity and the magnetic field strength, 
which supports the observations of Christensen et al. (2009). 
Future simulations of dynamos with density stratification and 
fully convective interiors will address the remaining critical fac- 
tors. 
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Appendix A. Simulation data 
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Table A. 2: Simulations results. The Prandtl number Pr is unity. The Pm values of bistable multipolar dynamos are marked with "*". 
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0.76 
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0.012 


0.25 






i 
i 


Ann 

OUU 


o 

/. jU 


70 
yj. 1 L 


U.Uuz 


u.uzzo 


U.40 


U.UUjV 


n o 
u.z 






i 
i 


7^n 

/ JU 


J. 1/ 


1 1 Q OQ 


U.Uo 1 


M H07^ 


U.o 1 


U.UIUJ 


u.z / 






U.J 


lUJU 


7 QQ 
j.yy 


oy. jj 


U. lUO 


U.U4J 1 


U.J J 


U.Uloo 


M 0/1 


1 x 10 


_4 


1 


500 


2.09 


62.39 


0.034 


0.0013 


0.63 
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1 


600 


2.42 


75.82 


0.042 


0.0017 


0.92 


0.8004 


0.43 






1* 


600 


2.28 


96.8 


0.0365 


0.0062 


0.33 


0.0086 


0.29 






1 


800 


2.99 


98.48 


0.058 


0.0023 


1.51 
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0.46 






1* 


800 


3.11 


130.97 


0.057 


0.0075 


0.63 


0.018 


0.32 






0.5 
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3.6 


64.81 


0.078 


0.0028 


0.80 


0.8891 


0.4 






0.5* 


1000 


3.74 


84.98 


0.066 


0.011 


0.43 


0.0111 


0.34 






0.5 
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109.98 


0.131 
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